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Abstract. The possible existence of strange stars in the universe will help in the understanding of various prop- 
erties of quantum chromodynamics, like asymptotic freedom and chiral symmetry restoration, which is otherwise 
very difficult to prove in laboratory experiments. Strange star properties were calculated using large color ap- 
proximation with built-in chiral symmetry restoration. A relativistic Hartree Fock calculation was performed 
using the Richardson potential as an interquark interaction. This potential has the asymptotic freedom and a 
confinement-deconfinement mechanism built into it and the present calculation employs an application of this 
potential with modified two scale parameters A and A', to find a new set of equations of state for strange quark 
matter. The linear confinement string tension from lattice calculations is 350 MeV and the Coulomb -like part has 
the parameter 100 MeV from deep inelastic scattering experiments. We also consider the effect of temperature, 
on gluon mass in a simple way, in addition to the usual density dependence. We obtained a set of new equation of 
states for strange quark matter and also found that the transition temperature from hadronic matter to strange 
matter is at 80 MeV, close to the 100 MeV estimated in litarature. Formation of strange stars may be the only 
signal for formation of quark-gluon plasma with asymptotic freedom and chiral symmetry restoration and this 
may be observable through many processes -such as for example through delayed 7 ray afterglow. 
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1. Introduction 

Two decades have passed since the existence of strange 
matter and strange stars was proposed by Witten (1984), 
and it is still difficult to prove or disprove the existence 
of a state of strange quark matter in its purest form. 
Besides the mass radius relation, efforts have been made 
to find other signatures for the existence of quark stars. 
The existence of absorption bands at 0.35 keV with har- 
monics at 0.7, 1.05, 1.4 etc, in IE 1207-5209 was in- 



terpreted in Sinha et ah (2003) as evidence for strange 
stars. It was then observed that six stars show harmonic 
band emission at energies corresponding to the same com- 
pressional modes at which the absorption bands were in- 
voked play et al., 2004j ). The SPI spectrometer on the 
INTEGRAL satellite detected a bright 511 keV line in the 
7 - ray spectrum from the bulge of the galaxy with a spher- 
ically symmetric distribution (Knodl eseder et al., 20 03). 
This has stimulated research to the fundamental physics 
that describes cosmological dark matter. Several authors 
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HBoehm et al., 2004] Oaknin and Zhitnitsky, 2005) have 
suggested that the line comes from positron annihilation, 
possibly in strange stars made of antiquarks. The impli- 
cation of cosmic separation of phases (Witten, 1984) is 
that strange stars (SS) can exist since the early universe 
and can themselves explain dark matter. A crucial con- 
dition for the survival of these stars, given by Alcock 
and Olinto (1989), is a large value of surface tension 
~ (178 MeV) 3 which can be obtained from strange star 
equation of state (EOS) (Bagchi et al., 2005). In his re- 
cent review, Weber (2005) argued that the constraints on 
mass radius relations from observational data support that 
the compact objects SAX J1808.4-3658, 4U 1728-34, 
RX J1856.5— 3754 and Her X— 1 are strange stars rather 
than neutron stars and on the basis of its low value of 
temperature, PSR 0205+6449 may be considered as a 
strange star. He also discussed the possibility of explaining 
some astrophysical phenomena like X-ray bursts, gamma- 
ray bursts and soft gamma repeaters with a strange star 
model. The most controversial is the isolated compact star 
RX J1856.5— 3754, which has a featureless thermal spec- 
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trum, and its X-ray data fits quite well with a blackbody 
distribution. With a distance estimate of 120 pc, Drake 
et al. (2002) gave the radiation radius of the star to be 
4.12±0.68 km. This indicates that the stellar radius is 4 
km to 8 km. Later, with a better estimation of the dis- 
tance of 175 pc, the estimated radiation radius R°° ~ 8 
km. However, there are controversies regarding the sur- 
face temperature of this star. With an estimated surface 
temperature (T°°) of 63. 5 eV, R°° ~ 4.4 km, and with 
T°° ~ 32 e y, R°° ~ 8 km flBurwitz et al., 2003| ). 

There is evidence, as recently discussed by Bombaci, 
Parenti and Vidana (2004) - for the existence of com- 
pact strange and neutron stars with larger radius and a 
shift from the latter to the former can account for delayed 
gamma ray bursts. 

PSR 1913+16 is the first pulsar to be discovered in a 
double neutron star system is, also known as the Hulse 
Taylor pulsar (1975). It is estimated to have a mass of 
1.441 Mq and the recently discovered pulsar (Lyne et al., 
2004) PSR J07373039A, forming a binary pulsar system, 
has a mass of 1.337 M . Thus the mass of pulsars may 
lie near these values, the observed variation of periodic- 
ity of the pulsars indicates the variation in their masses. 
Although it may be a stringent lower-limit constraint for 
neutron stars masses, which are bound mainly by gravi- 
tational forces, the same does not hold for compact ob- 
jects, which are self bound, like a system of quark matter. 
Theoretically there is no physical reason to have a lower 
limit for the mass of strange stars, where strong interac- 
tions play an important role in stellar binding. The same 
EOS will allow even a very small quark nugget. 

In the literature, there are several EOSs for 
strange quark matter, starting from the Bag 
model dAlcock et al., 1986| |Haensel et al, 1986| 
|Kettner et al., 1995| ) to the recent models like the mean 
field model with interacting quarks (De y et al., 19 98), 
the strange star EOS model in the perturba- 
tive quantum chromodynamics (QCD) approach 
< |Fraga et al., 20011, those of the chiral chromodi- 
electric model JMalheiro et al., 2003| ) or Dyson-Schwinger 
model dBlaschke"et"^L7"l999| l, etc. In our present work, 
we will consider only the relativistic mean field model by 
Dey et al. (1998). 

In the relativistic Hartree-Fock calculation of a 
strange quark matter EOS to form a compact star 
(De y et al., 1998| , a phenomenological inter-quark in- 
teraction, namely the Richardson potential, was used. 
Strange stars are more compact than conventional neu- 
tron stars and fit into the Bodmer (1971) - Witten (1984) 
hypothesis for the existence of strange quark matter. They 
were able to produce a very compact stellar structure with 
a mass of about 1.44 Mq and radius of about 7 km. The 
original Richardson (1979) potential was designed to ob- 
tain the mass spectrum of heavy mesons (Charmonium 
and Upsilon). It includes two features of the qq force, the 
asymptotic freedom (AF) and confinement, however, with 
the same scale, A. 



The A parameter needed was 400 MeV. It was further 
applied to light meson spectroscopy and baryon proper- 
ties with the same value of A (Crater and van Alstine, 
1984, Dey, Dey and Le Tourneux, 1986). Although both 
the confinement and AF scales are chosen to be the same 
in the Richardson potential, it is the confinement, not the 
AF part, that plays the more important role in determin- 
ing the static properties of hadrons. However, when the 
AF part is important (as in strange star property calcu- 
lations), A would be much smaller, around 100 MeV, as 
it is the scale for AF, obtained from perturbative QCD. 
Indeed, for self bound high density strange quark matter, 
Dey et al. (1998) found that A needed to be ~ 100 MeV. 
In the present paper, proper scales are chosen for both the 
parameters, as determined from a calculation by Bagchi 
et al. (2004) by fitting the masses and recently determined 
magnetic moments of fl~ and A ++ . 

In section we describe the basic formalism of the 
model as given in Dey et al. (1998) and the calculations 
incorporating finite temperature effects. In section [31 we 
present a fine tuning of the Dey et al. (1998) model, with 
a modified version of the Richardson potential, and also 
incorporate a temperature dependence of gluon mass. In 
section 0] we show the results and the plots, and draw 
conclusions and summarize in section [S] 

2. The model 

There are no free parameters in QCD with which one 
can systematically build up a perturbation series. So 



't Hooft (1974) suggested use of the inverse of the number 
of colors - which is the key quantum number in QCD - as 
an expansion parameter. In this case a baryonic system 
like stars can be explored in relativistic tree level calcu- 
lations with a potential that may be obtained from the 
meson sector phenomenology (Witten, 1979) In this con- 
nection the Richardson potential is a very good choice for 
interquark interactions. 

The original Richardson potential is given by: 
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At high density (high Fermi momentum) quarks are 
deconfined as the bare potential is screened. Confinement 
is softened and the AF part takes over. This bare potential 
in a medium will be screened due to gluon propagation and 
(kj-kj) 2 of equation Q is replaced by [(kj — kj) 2 +£>~ 2 ]. 
The inverse screening length, which can be identified 
with the gluon mass to s , to the lowest order is : 



{D -x ? = 2ao 
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where k\ , the Fermi momentum of the i-th quark is ob- 
tained from the corresponding number density: 

k{ = Ktt 2 ) 1 / 3 (3) 

and ao is the perturbative quark gluon coupling. 
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In this model, chiral symmetry restoration at high den- 
sity is incorporated by assuming that the quark masses are 
density dependent : 



Mi = rrii + MQsech 



( n B 



\Nn 



i = u,d,s 



(4) 



where tib = {n u +nd+n 3 )/3 is the baryon number density; 
n = 0.17 fm is the normal nuclear matter density; n u , 
rid, n s are the number densities of u, d and s quarks re- 
spectively and N is a parameter. The current quark masses 
(rrii) are taken as : m u = 4 MeV, = 7 MeV, m s = 
150 MeV . It is ensured that in strange matter, the chem- 
ical potentials of the quarks satisfy (3 equilibrium and 
charge neutrality conditions. The parameters Mq and iV 
are adjusted in such a way that the minimum value of E/A 
for u,d,s quark matter is less than that of Fe 56 , so that 
u,d,s quark matter can constitute stable stars. The mini- 
mum value of E/A is obtained at the star surface where 
the pressure is zero. The surface is sharp since strong in- 
teraction dictates the deconfinement point. An electron 
cloud forms above this sharp surface - since they are not 
held by strong interaction. However, the minimum value 
of E/A for u-d quark matter is greater than that of Fe 56 
so that Fe 56 remains the most stable element in the non- 
strange world. With the obtained EOSs, they solved the 
hydrostatic equilibrium equations (Tolman-Oppenheimer- 
Volkov equation) to get the structure of the stars. 

Finite temperature T was incorporated in the system 
through the Fermi function (Ray et al., 2000b): 



FM(k,T) = 



1 



e (e-e F )/T + 1 

with the flavour dependent single particle energy 



y/k 2 +Mi(p) 2 + Ui(k). 
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where g is the spin-colour degeneracy, equal to 6 and I is 
the number density for <f)(e) = 1 and energy density for 
<t>{e)=e. 

The system is highly degenerate even at very high T 
which is around 70 MeV since the chemical potential is 
very high, of the order of several hundred MeV . The en- 
tropy was calculated as follows: 

o poo 

S{T) = - / k 2 [FM(k, T)ln(FM(k, T)) 

JO 

+ (1 - FM(k, T))ln(l - FM(k, T))]dk (8) 
The pressure was calculated from the free energy 
/ = e - Ts 
as follows: 

fif- 

(9) 



^ dpi 



Ray et al. (2000b), found that stable star structure 
cannot be obtained with a temperature greater than 70 
MeV. But here the effect of temperature on gluon mass 
was not considered. Also, the parameter A needs to be 
different. Can one separate the confinement and AF part 
and apply the same two parameters to hadron properties 
as well as to high density quark matter? Moreover, what 
happens to the strange star properties if the effect of tem- 
perature on gluon mass is considered? 

Successful attempts have been made to construct a 
two parameter Richardson potential to explain energies 
and magnetic moments of baryons. Recently the magnetic 
moments of A ++ and Or with three u and three s valence 
quarks respectively have been accurately found in exper- 
iments. These sensitive properties are fitted with a modi- 
fied Richardson potential with different scales for confine- 
ment (~ 350 MeV) and AF (~ 100 MeV) (Bagchi et al., 
2004) . It was also shown that baryonic properties depend 
much on the confinement part and less on the AF part as 
baryons are confined quark systems. It is natural to apply 
this new potential to star calculations since they are con- 
strained by baryonic data. Retaining the approach of Dey 
et al. (1998), we used this modified Richardson potential 
to derive our EOS. 

3. Modifications 

The properties of AF and confinement of quarks in the 
model is embedded in the potential we have used, i.e., the 
Richardson potential. We have extensively searched for 
the two parameter form of this potential that would sepa- 
rately reflect the properties of AF and quark confinement. 
In the light of the new modifications, the Richardson po- 
tential takes the form: 
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with A' the confinement property and A asymptotic free- 
dom. 

The term ( q2;„ (1 +q2 /a2) - 3* J is asymptotically zero 



v Q 2 in(l+Q 2 /A 2 ) 

for large momentum transfer Q 2 

A'' 



and the 



term ^ explains the confinement reducing to a linear 
confinement for small Q 2 . 

Previously Bagchi et al. (2004) used the same form of 
the potential to fit baryonic properties and evaluated the 
appropriate values of A and A'. We have borrowed the 
same values (A — 100 MeV and A' ~ 350 MeV) in our 
present calculation. 

Ray et al. (2000b) calculated finite temperature effects 
on strange matter and strange stars, using the original 
one parameter Richardson potential. There, the temper- 
ature dependence of gluon mass was also not considered. 
We have incorporated this effect following Alexanian and 
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Nair (1995), who obtained a magnetic screening for the 
quark-gluon plasma at high density after performing self- 
consistence summation of one loop self-energy. It has the 
form M w (2.384)^1 w here C = N for SU(N) gauge 

2 

theory and g is the strong coupling constant |— = a®. 
Adding this term to the usual density dependent one, the 
inverse screening length becomes: 

(£) - 1)2= 2^o j2 k^iklY+m* + 7.152 a T (11) 

i— u,d,s, 

where we have used our standard notations as before. We 
have incorporated this form of the gluon mass (inverse 
screening length) for calculations of our new sets of SS 
EOSs. 

4. Results 

We check our new potential by calculating strange star 
properties at zero temperatures. We obtain different EOSs 
for different set of parameters tabulated in Table d For 
comparison, we also present one EOS SSI obtained by Dey 
et al. (1998) using the unmodified potential in Tabled In 
Figured we compared the M-R plots for three of our new 
EOSs with SSI. The compactness (M/R ratio) of the stars 
has increased for the new EOSs than compared to SSI. 

From Tables d and the nature of EOSs and M - R 
relations do not change appreciably when using the modi- 
fied potential if we adjust other model parameters suitably 
(see Figured- 

In our new EOSs, we have kept the value of the param- 
eter N to be the same (except EOS C, where we changed 
N to see the effects). However, the value of the param- 
eter «o for our new EOSs is significantly different from 
SSI. Here, we chose a larger value because, as the confine- 
ment parameter has now increased from the previous value 
of 100 MeV to 350 MeV, ao should be increased to in- 
corporate an earlier onset of deconfinement. The present 
value of a is 0.55-0.65 compared to 0.20 of Dey et al., 
(1998). Such a value is in agreement with theoretical re- 
sults ( |Ray et al., 2000a| . From Tabled w e infer that the 
parameter TV has the greatest effect on the star structure 
and Mq the least. 

With the parameters used to obtain EOS F, we have 
calculated strange star properties at finite temperatures. 
We have found that stable stellar structure can be ob- 
tained up to a temperature of 80 M eV which is closer to 
Witten's conjectured temperature of cosmic separation of 
phases (100 MeV). Our results are shown in Figuresdand 
d The figures imply that in our strange star, the density 
(Pb) cannot be less than 4 po because negative pressure is 
physically not feasible. This means that the strange star 
has a sharp surface at a density ps ~ 4 po where the pres- 
sure goes to zero. This property of the presence of sharp 
surfaces in strange stars has also been seen for other mod- 
els, like the bag model, the chromo-dielectric model, etc. 

In Figure 01 we have shown both the density and tem- 
perature dependence of gluon mass. It is clear that at a 
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Fig. 1. Mass - Radius relations for different parameter 
sets. Dashed curve is for EOS C, dotted curve is for EOS 
B, solid curve is for EOS SSI and the dot - dashed curve 
is for EOS F. 
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Fig. 2. EOSs at different temperatures, all with the same 
parameters as EOS F. The first vertical line from the left 
shows the location of the zero pressure (and hence stellar 
surface) for the EOS at 80 MeV , which is around 3.45. 
The second vertical line is for an EOS at 30 MeV, where 
the zero of the pressure is at 5.05. 

given temperature, gluon mass increases with density and 
at a given density, gluon mass increases with temperature. 

5. Conclusions and summary 

A set of new EOSs for strange matter is presented - both 
for zero and very high T, using the Richardson potential 
with the value of A' ~ 300 to 350 MeV and the value 
of A = 100 MeV, the two scales for confinement and 
AF respectively. This explains both the properties of the 
deconfined quark matter (our present work) and the prop- 
erties of confined 3u and 3s baryons (Bagchi et al., 2004). 
Although the confinement part A' is stronger, leading 
to a sharp surface, it is softened by medium effects, de- 
veloping a screening length. Inside the star, the AF part 
is more important. We have found that for a wide range 
of parametric variation, the strange matter EOS gives a 
minimum of energy per baryon (E/A) that is lower than 
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Table 1. Data for different EOSs, each having two different values for A and A'; A is always taken as 100 MeV. 
Change in the confinement scale A' does not change the minimum energy per baryon significantly. 



EOS 


A' 


M q 


N 


Q() 


(E/A) mm 




R 


Label 


(MeV) 


[MeV) 






[MeV) 




(km) 


A 


350 


325 


3.0 


.55 


874.0 


1.532 


7 All 


B 


350 


325 


3.0 


.65 


907.0 


1.470 


7.140 


C 


300 


325 


2.7 


.55 


877.0 


1.554 


7.530 


D 


300 


325 


3.0 


.55 


906.0 


1.463 


7.110 


E 


300 


335 


3.0 


.55 


912.0 


1.447 


7.022 


F 


350 


345 


3.0 


.65 


920.3 


1.436 


6.974 



Table 2. One EOS (SSI) obtained by Dey et al.. (1998) using an unmodified Richardson Potential with both A and 
A' equal to 100 MeV. 



EOS A M q N a (E/A) mm M max R 

Label (MeV) (MeV) (MeV) (Mq) (km) 

SSI 100 310 3 .20 889.0 1.437 7.06 




(E/A) Fe56 = 930.6 MeV which ensures that the system 
is stable, and unlike neutron star - like structures, our 
strange stars are not solely bound by gravitational forces 
but also by strong forces. Moreover, considering tempera- 
ture dependent screening of the potential, we have found 
that strange stars can sustain stable configurations up to 
a temperature of 80 MeV. 

The indirect proof of the existence of strange stars 
is becoming more prominent because there are many ob- 
served features from compact objects that cannot be ex- 
plained by the standard neutron star model, like spectral 
features and the controversies built around the star RX 
J1856.5— 3754. Another feature, superbursts, faces diffi- 
culties in using the standard carbon burning scenario in 
the neutron star crust. The carbon, nitrogen and oxygen 
mass fraction (Zqno) must be larger than Zqno,q, where 
the latter refers to the standard value found in the sun 
(Cooper et al. 2005). Sinha et al. (2002) predicted that an 
alternative to the carbon burning scenario in a neutron 



star to explain superbursts, can be formation of diquark 
pairs on the strange star surface and their subsequent 
breaking, thus giving a continuous and prolonged emis- 
sion of energy that can be comparable to the superbursts. 
Recently, Page and Cumming (2005) tried to explain su- 
perbursts as carbon burning of matter accumulated on the 
electron cloud of a strange star. There are many other fea- 
tures that can be fitted better with a strange star than a 
neutron star model. So, the study of strange matter and 
strange stars is important and the formulation of a realis- 
tic strange star model is needed. 

Acknowledgements. MB, MD and JD acknowledge DST grant 
no. SP/S2/K-03/2001, Govt, of India and the IUCAA, Pune, 
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